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Abstract. Using the filled bolometer array SHARC-II on the Caltech Submillimeter Observatory (CSO), we have obtained 
deep (arms ~ 15mJybeam^^), high-quality 350-/im maps of five of the most luminous high-z radio galaxies known. In all 
cases the central radio galaxy is detected at the > 3a level, and in some cases the high resolution of SHARC-II (FWHM = 9") 
allows us to confirm the spatially extended nature of the dust emission. In PKS 1138—262 (z = 2.156), 8C I909-I-722 
(z = 3.538) and 4C41.17 (z — 3.792), additional sources - first discovered by SCUBA at 850-/im and believed to be dusty, 
merging systems associated with the central radio galaxy - are detected at 350-^m. Eurthermore, in PKS 1138 and 4C 41. 17 
additional SHARC-II sources are seen which were not detected at 850-/im, although the reality of these sources will have to be 
confirmed by independent submm observations. Thus, our observations seem to support the notion of extended star formation 
taking place in radio galaxies at high redshifts, and that these systems sit at the centers of overdense regions harbouring a 
complexity of dusty and vigorously star forming systems. At the redshift of the radio galaxies, the 350-/im observations 
sample very close to the rest-frame dust peak (typically at ~ 100- ^m), and they therefore contribute a particularly important 
point to the spectral energy distributions of these sources, which we use this in conjunction with existing (sub)millimeter data 
to derive EIR luminosities, dust temperatures and spectral emissivities of the central radio galaxies. 
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1. Introduction 

High-0 Radio Galaxies (HzRGs) are typically identified 
in low frequency radio surveys as Ultra Steep Spectrum 
(USS) radio sources (De Breuck et al. 2000), which in con- 
junction with near-IR colour criteria provide a highly ef- 
ficient selection scheme. The most powerful [PnsAiHz 
10^® /i^^ W Hz^^) HzRGs have become an invaluable tool 
for testing models of galaxy formation, as radio galaxies are 
amongst the most massive galaxies at all redshifts, and as 
such they trace out the most biased peak in the primordial 
density field. Furthermore, owing to their accurate radio po- 
sitions and morphologies, as well as their extreme properties, 
such as large (optical) luminosities, extended X-ray and Lya- 
halos, HzRGs are relatively easy to study compared to other 
populations of high-z galaxies. 

The first pointed submillimeter (submm) observations of 
high-z radio galaxies demonstrated that a large fraction of 
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such systems are extremely FIR luminous (Lpm ^ 10^^ Lq ) 
and contain large amounts of dust (Md ~ 10^ A/q - Dun- 
lop et al. 1994; Chini & Krugel 1994; Ivison et al. 1995). 
Subsequently, systematic 850-/im surveys with SCUBA of 
HzRGs by Archibald et al. (2001) and Reuland et al. (2004) 
confirmed these initial findings, and showed that the dust con- 
tent of HzRGs increases strongly with redshift. The latter was 
interpreted as HzRGs being increasingly younger (younger 
stellar populations) at higher redshifts, since while the FIR 
luminosities can be powered by either Active Galactic Nu- 
clei (AGN) or starburst activity, the large dust masses re- 
quires that vigorous star formation has recently been taking 
place. Thus, even in HzRGs with their powerful AGNs, star- 
burst activity is likely to contribute substantially to the ex- 
treme FIR luminosities, albeit probably less so than in high- 
z submm-selected galaxies which are predominantly star- 
burst dominated (Alexander et al. 2003). Furthermore, about 
a handful of HzRGs are known to harbour large reservoirs 
of molecular gas (M(H2) ^ 10^^ Mq - e.g. Papadopoulos 
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et al. 2000; De Breuck et al. 2003; Greve et al. 2004; De 
Breuck et al. 2005; Klamer et al. 2005) - enough to fuel a 
~ 1000 Mq yr^i starburst for ~ 10* yr. Finally, the peak of 
the dust and CO emission is often not centered on the central 
radio galaxy, but offset and spatially extended on scale of up 
to 30-200 kpc. Some of the most striking and visually grati- 
fying examples of this came from submm (850-/xm) SCUBA 
observations of seven high-z radio galaxies, which not only 
revealed extended 30 kpc) dust emission from the central 
radio galaxy itself, but also showed an overdensity by a factor 
of two compared to the field of sub-mm sources around the 
radio galaxies (Stevens et. al 2003), suggestive of extended 
star formation on hundreds of kilo-parsec scales. 

2. Results 

Here we present 350-/im Ussajou-style scan-maps obtained 
with SHARC-ii (Dowell et al. 2002) on the Caltech Submil- 
limeter Telescope (CSO) on Hawaii of 5 of the 7 HzRGs 
mapped with SCUBA at 850-/xm by Stevens et al. (2003). 
The motivation for doing this was two-fold: 1) to utilise the 
superior resolution of SHARC-ii (fwhm = 9") over that of 
SCUBA (FWHM = 15") to confirm the spatially extended dust 
emission from the HzRGs, and the reality of the many com- 
panion objects reported by Stevens et al. (2003), and 2) to 
determine accurate FIR luminosities and dust masses in these 
objects. SHARC-II observations are particular useful for the 
latter since at the typical redshift of the HzRGs (z ~ 3.5) they 
probe the ~ 100 fim rest-frame peak of the dust emission. 
Stevens et al. targeted these 7 sources because they were the 
brightest sources at 850-/im of the SCUBA observed sample 
of HzRGs by Archibald et al. (2001). Thus, it is important to 
keep in mind that our samples represents the most luminous 
and extreme examples of radio galaxies. 

PKS 1138-262 At 350-/Ltm the central radio galaxy (ab- 
breviated RG, and marked by a cross in Figure 1) is seen to 
be associated with a ~ 5cr emission feature consisting of two 
connected knots roughly aligned with the axis defined by the 
kiloparsec scale radio jets (marked with tick marks around 
the central RGs in Figure 1). This double-knotted structure 
is most hkely real since it is also discernable at 850-/xm, 
and in both maps the south-western most source {S^^Ofim = 
(92 ± 18) mJy, 7" south-west from the RG) is the bright- 
est, see Table 1. However, at the position of the RG, which 
coincided with the weaker knot (S'350;um = (61 ± 15) mJy), a 
1200-/im flux density of (1.5 ± 0.5) mJy was detected using 
MAMBO on the IRAM 30 m Telescope (Carilh et al. 2001). 

The brightest SCUBA source (SgsoMm = (7-8 ± 2.2) mJy) 
detected by Stevens et al. in this field lies about 70"east of 
the RG and coincides with a 1.9cr peak in the 350-/Ltm map. 
Several explanations may account for the lack of a significant 
350-/im detection of this very bright SCUBA source: a) the 
850-/xm detection could be spurious (it lies at the edge of the 
SCUBA map), b) the source lies near the edge of the sharc- 
II map where the noise is higher which makes it less likely 
to be detectable, or c) the source is real, but at such a high 
redshift that the 350-/im band is shifted beyond the dust peak 



and samples the SED so far down the Wien tail that it is un- 
detectable given the depth of the maps. We are able to put a 
3(7 upper limit on its 350-/im flux of Ss50iJ.m ^ 95 mJy. 

Another, much weaker SCUBA source, ^ 37" west of the 
RG, coincide with ~ 2a peaks in the 350-/im map. We de- 
rive an upper hmit on its 350- ^m flux of S'sso^m ^ 68 mJy. 
Finally, a third SCUBA source ^ 57" south of the RG falls 
outside the 350-/im map. Thus, while we verify the double- 
knot structure of the central radio galaxy, we do not confirm 
any of the apparently associated sources seen with SCUBA. 

8C 1909+722 The 350-iJ,m sharc-ii map shows a strong 
(^ 5cr) central source (Sssoum = 90 ± 15mJy), and a 
second significant (> 4(j) source about 38"south-eastwards 
iS350fj,m ~ 148 ± 32mJy). The former is offset from the 
RG-position by ~ 6", while the latter lies within ~ 8" of a 
strong SCUBA source detected about 46"south-east of the RG 
(Stevens et al. 2003). These offsets can easily be accounted 
for considering that a) we have not corrected for any poten- 
tial astrometrical offsets between the SHARC-II and SCUBA 
maps, 2) the beam widths are large (fwhm ~ 15" in the case 
of SCUBA), and that c) the two sources are clearly extended 
at 850-^m. 

The RG is marginally resolved at 350-/im, and exhibit the 
same general shape as seen at 850-/im, thus apparently con- 
firming the claim by Stevens et al. (2003) that the (850-/im) 
dust emission in this object (and its brightest companion) is 
extended on tens of kpc scales. This seems to be backed up 
by the detection of a weak (^ 4(t, S'sso/^m = 70 ± 17mJy) 
source only ^ 19" north-east of the RG, and thus within the 
north-east extended 850-/im emission from the RG. 

A third much fainter SCUBA source, ~ 35" west of the 
RG and within the SHARC-II map, is not detected at 350- /[xm, 
and an upper flux limit of Ss^o/^m ^ 28 mJy is derived. 

4C 60.07 This source received by far the least amount of 
integration time of the 5 HzRGs observed with SHARC-II, and 
as a result the source is only just detected at the ~ 3a level 
(S350nm = 59±22 mJy). This source was one of the strongest 
sources detected at 850-^m by Stevens et al. and even showed 
some, albeit marginal, evidence of being extended. A faint 
850-/im source found ^ 30" north of the RG coincides with 
extended 'da emission in the SHARC-II map. However, given 
its small 850-/um flux (S'aso^m = 5.9 ± 1.3 mJy), and the 
somewhat larger noise level in that part of the SHARC-II map, 
we do not consider this source robustly detected at 350-/:im. 
We derive a 3a upper flux limit of ^sso^m ^91 mJy 

4C 41.17 With 5 sources detected at > 5a significance, 
this is by far the richest field observed at 350-/xm. All three 
SCUBA sources detected by Stevens et al., including the RG 
which even appears to be slightly extended, are detected at 
350-/im, and their flux densities are listed in Table 1. Ivison 
et al. (2000) reported an additional three faint SCUBA sources. 
However, only one of those sources lies within the SHARC-II 
map and is not detected to a 3a limit of S'sso^m ^ 67 mJy. 
The most significant (~ 8a) 350-/um source lies ~ 53" 



3 SPECTRAL ENERGY DISTRIBUTIONS, FIR LUMINOSITIES AND DUST MASSES 
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Fig. 1. SHARC-II 350-/im signal-to-noise maps of the five radio 
galaxies observed. Contours are at 3, 4, 5, 6a, where a is the rms 
noise in each map, see the individual panels. The SHARC-II beam at 
350-/im (FWHM — 9") is shown in the lower right corner of each 
panel. The crosses mark the position of the radio core. The open cir- 
cles and the associated labels indicate the 850-^tm sources detected 
with SCUBA by Stevens et al. (2003) and Ivison et al. (2000). The 
tick marks around the central radio galaxy indicate the alignment of 
the large-scale radio jets, and the large circles outline the extend of 
the SCUBA maps of Stevens et al. (2003). 
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Fig. 2. Radio-FIR/submm spectral energy distributions of the 5 
HzRGs observed with SHARC-II. The dust temperature and spec- 
tral index were allowed to vary freely in the fitting process, except 
in the case of 8C 1909+722 where /? was fixed at 1.0 due to the 
small number of data point. The radio part of the spectrum was fitted 
with a parabola (dashed curves). The combined radio to FIR/submm 
SEDs are shown as solid curves. 

8C 1435+635 Despite being the faintest radio galaxy at 
850-/xm, the central radio galaxy is clearly detected (> 5cr) 
in the SHARC-II maps (Ss^o^m = 51 ± 9mJy). Similar to 
Stevens et al., who failed to resolve the galaxy at 850-/im, 
there is no evidence of extended emission at 350-/im either, 
suggesting that this source may be relatively compact. The 
only other SCUBA source detected in the 8C 1435+635 field 
lies about 60" north-west of the RG and coincides with faint 
2(7 emission at 350-/im. We derive an upper flux limit of 
S350iJ.m H 68mJy. 



south-west of the RG, and although it is not detected at 850- 
/im it appears to be associated with the very bright SCUBA 
(and SHARC-II) source ~ 30" north-east of it. The two 350- 
fim sources malce up a large, extended structure aligned along 
the jet-direction. The reality of this structure, and of the most 
south-western source in particular, will have to be confirmed 
by independent submm observations. Finally, a ^ 5cr source 
is detected at 350-/im, about 68"north of the RG. Unfortu- 
nately, the source lies outside the SCUBA map and thus can- 
not be confirmed at 850-/xm. Given that this source lies right 
at the map edge where the noise is higher, we suspect that it 
may be spurious. 

Although, we cannot exclude that some of the SHARC- 
II sources (in particular the ones closest to the map edges) 
are spurious, the confirmation of the three SCUBA sources 
reported by Stevens et al. (2003), supports their conclu- 
sions that the environment of 4C 4 1.1 7 has an overdensity of 
submm sources. 



3. Spectral Energy Distributions, FIR 
luminosities and Dust Masses 

For each of the central radio galaxies, an optically thin grey- 
body spectrum of the form 5*1, oc jy3+/3('g'ii>/*:Td _ 1)~^, was 
fitted using all FIR/submm data points available in the liter- 
ature. In all cases but 8C 1909+722 more than one indepen- 
dent 850-/im SCUBA measurement existed (Ivison et al. 1995, 
1998; Hughes et al. 1997; Archibald et al. 2001, Reuland et 
al. 2004), all of which were included in the fitting-process, 
and in a few cases even 1200-, 450- and/or 350-/im observa- 
tions had been made. In the case of 8C 1909+722 where the 
FIR/submm SED has been sampled at only two separate fre- 
quencies we allowed only Td to vary in the fit for fixed values 
of /3 = 1.0, 1.5 and 2.0. The best fit was found for /3 = 1.0. 
The resulting dust SEDs are shown as dotted curves in Figure 
121 and the fitted Td and (3 values are listed in Table 1 . The ra- 
dio part of the spectrum was fitted with a parabolic function 
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Table 1. Observed 350-, 450-, 850- and 1200-/xm fluxes and derived properties of the 5 HzRGs and their associated submm 
sources (which fafl within the SHARC-II maps). The numbering foUows that of Stevens et al. (2003) and Ivison et al. (2000). 



Source 


z 




cfc 




'-'1200/J,m 




/9 




Md 






[mJy] 


[mJy] 


[mJy] 


[mJy] 


[K] 




[xlO'^'L©] 


[xIO^Mq] 


PKS 1138-262a 


2.156 


61 + 15 




3.7 + 2.4 


1.5 + 0.5"= 


49 + 6 


1.9 + 0.2 


2.1+0.5 


2 + 1 


PKS 1138-262b 


2.156 


92 ± 18 




3.8 + 2.4 












1. 




< 95 




7.8 + 2.2 












3. 




< 68 




2.2 + 1.4 












8C 1909+722 


3.536 


90 + 15 




34.9 + 3.0 




69 + 3 


l.O'^ 


5.7 + 0.9 


13 + 2 


1. 




148 + 32 




23.0 + 2.5 












3. 




< 28 




4.3 + 2.1 












4. 




70 + 17 
















4C 60.07 


3.788 


59 + 22 


69 + 23 


23.8 + 3.5 




58 + 8 


1.0 + 0.2 


2.9 + 1.1 


15 + 6 


1. 




< 91 




6.3 + 2.1 












4C41.17 


3.791 


71 + 17 


35.3 + 9.3 


12.0 + 2.3 


2.7 + 0.3 


60 + 3 


1.3 + 0.1 


4.1+0.9 


8 + 2 


1. 




108 + 17 


34.1 + 9.3 


21.2 + 2.9 












2. 




184 + 34 




6.2 + 1.9 












4. 

6. 




< 67 
197 ± 24 


< 26 


2.8 + 0.8 
< 5 












7. 




181 + 40 
















8C 1435+635 


4.261 


51 + 9 


23.6 + 6.4 


6.0 + 2.1 




64 + 5 


1.4 + 0.1 


3.8 + 0.7 


4+1 


1. 




< 68 




3.9 + 1.8 













°- The total iltix density wittiin a r — 15" aperture. The flux errors include the ~ 15 per-eent calibration uncertainty; 450-^ni flux densities are from Archibald et al. (2002) and Ivison et al. (2000); ^ 850-;xni flux densities are front Stevens 
et al. (2003), except for PKS 1138 — 262a and b which are our own measurements, and 4C41. 17-4 which is from Ivison et al. (2000); 1200-^m flux densities are from Chini & Krugel (1994) and CarilU et al. (2001); ^ Fixed value. 



(dashed curves), and the combined FIR/submm-radio SED is 
shown as solid curves in Figure|2] 

Due to the fact that our 350-/im observations probe the 
SEDs near their peaks at ~ lGO-/im rest-frame wavelengths, 
we are able to constrain dust temperatures and FIR lumi- 
nosities much better than previously. In general, the derived 
dust temperatures are rather high (lowest value is = 
49 K), compared to the median dust temperatures of submm- 
selected galaxies (Td ^ 36 K- Chapman et al. 2004). While 
HzRGs, in contrast to SMGs where the AGN rarely plays 
a dominant role (Alexander et al. 2003), harbour powerful 
AGN, capable of heating the dust in the nuclear regions to 
high temperatures, it is difficult to envisage them doing so 
on scales of > lOkpc as suggested by the extended 350- 
/Ltm emission. The spectral emissivity index varies between 
/3 1 — 2, a similar range to what is observed in SMGs. We 
note that for objects at z ~ 3, /3 is much better constrained by 
850- and/or 1200-/im observations which probe further down 
the Rayleigh-Jeans tail than 350-/im. 

FIR-luminosities obtained by integrating the FIR/submm 
SED over the wavelengths range 200 — 1000 /im range are 
listed in Table 1, along with dust masses derived from the ob- 
served 350- /j,m fluxes using a mass absorption coefficient of 
K(yr) = G.ll(i'r/352GHz)'^ in units of m^kg^^. Assuming 
that only half of the FIR-luminosity in these objects is due 
to star formation, we find that on average they form stars at 
a rate of 1500 A/q yr^^. The large dust masses suggest 
that star formation must have been going on for some time 
in these systems. That such large star formation rates can be 
sustained for a long enough time to form a giant elliptical 
(-- IGyr) is shown in the cases of 4C 60.07 and 4C41.17 
where the detection of CO has revealed large amounts of 



molecular gas (~ 10^^ Mq - Papadopoulos et al. 2000; 
Greve et al. 2004; De Breuck et al. 2005). 

4. Concluding Remarks 

To conclude, our SHARC-II maps confirm the over-densities 
of dusty galaxies detected in HzRG fields by SCUBA. Fur- 
thermore, the superior spatial resolution of these data allow us 
to confirm that dust emission associated with both the HzRGs 
and companions is extended on galaxy-wide scales. Our re- 
sults suggest that mergers play a key role in the formation of 
these galaxies. 
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